Penetration depth for shallow impact cratering 
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We present data for the penetration of a variety of spheres, dropped from rest, into a loose non- 
cohesive granular medium. We improve upon earlier work [Uehara et al. Phys. Rev. Lett. 90, 
194301 (2003)] in three regards. First, we explore the behavior vs sphere diameter and density more 
systematically, by holding one of these parameters constant while varying the other. Second, we 
prepare the granular medium more reproducibly and, third, we measure the penetration depth more 
accurately. The new data support the previous conclusion that the penetration depth is proportional 
to the 1/2 power of sphere density, the 2/3 power of sphere diameter, and the 1/3 power of total 
drop distance. 

PACS numbers: 45.70.Ht, 45.70.Cc, 83.80.Fg, 89.75. Da 
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The mechanics of granular media continue to defy our 
intuition. In spite of their ubiquity in everyday life and 
industry, we have no fully reliable rules for predicting re- 
sponse to an applied force 0,0,13. If the forcing is weak, 
the medium remains at rest and the local disorder in the 
packing gives rise to ramified force chains with structures 
much larger than the grain size. If the forcing is strong, 
the medium can flow. But when will the medium yield? 
Will the flow be smooth or intermittent? How do veloc- 
ity and density vary with position and time? No experi- 
mental characterization currently exists that can be used 
to predict response in all other sample and forcing ge- 
ometries. For example, flow down an incline offers little 
insight as to how the same medium would flow between 
rotating cylinders or on a vibrated plate. 

Recently we investigated the mechanics of impact by 
projectiles dropped into granular media 0, Q . This is 
a situation of natural interest. Some of us would like to 
understand how far our feet sink into the sand when walk- 
ing, running, or jumping at the beach. Others of us would 
like to understand the lie of our golf ball in a sand trap. 
Still others would like to know how far a warhead can 
penetrate the earth prior to detonation . It's also inter- 
esting to consider the effect of impact on the medium it- 
self: the nature of the granular splash Bja and the mor- 
phology of the resulting crater [ij [ToL 
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Our motivation is more general: to find a non-contrived 
situation permitting the unusual nature of granular me- 
chanics to be both highlighted and characterized. Projec- 
tile impact is ideal on both counts. It's interesting that 
penetration is nonzero even for near-zero impact speed, 
but grows only very slowly with projectile energy. The 
average stopping force (F) can be very small, but can 
also increase dramatically for stronger impacts. This un- 
usual mechanics can be studied quantitatively from the 
penetration depth d via a simple statement of energy 
conservation: 



(F) = mgH/d, 



(f) 



free-fall height h and the penetration depth d (see inset 
of Fig.QJ. 

In Ref. Q we measured the penetration of spherical 
projectiles of various densities p b , and diameters -Dt,, into 
loose non-cohesive granular media of various densities p g 
and angles of repose tan -1 p. In all cases the minimum 
free fall height was nearly zero, and the maximum pen- 
etration depth was comparable to the ball diameter. All 
of our data collapsed according to an empirical scaling 
relation, 
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where m is the projectile mass, g = 9.8 m/s 2 , and H is 
the total drop distance. Note that H is the sum of the 



In Ref. we showed how this naturally generalizes to 
cylindrical projectiles, independent of the tip shape. If 
true, Eq. J5J) has several interesting implications. First, 
it implies via Eq. JQ) that the average stopping force is 
proportional to the tangent of the repose angle, /i, con- 
sistent with the notion that it represents a friction coef- 
ficient 0. Second, it implies that the granular medium 
can be extremely fragile, suffering a nonzero penetration 
even for zero free-fall height. Eq. J2J) gives this minimum 
penetration as 

d a = (O.U/p)V 2 (p b / Pg f/ 4 D b . (3) 

The penetration depth formula then can be recast di- 
mensionlessly as d/d a = (H/do) 1 ^ 3 . Third, since Eq. (J2J 
is dimensionally complete, it suggests that the effects of 
grain-grain cohesion and interstitial air are both negligi- 
ble. If they were not, then even further physics would 
have to enter to cancel the extra units. Air and cohesion 
effects can also be ruled out because we found identical 
penetrations for granular media that are identical except 
for particle sizes [4J- This is to be expected, according to 
the Geldart classification scheme of fluidization behavior 
vs particle size and density (16)] . Fourth, and perhaps 
most curious, Eq. J2J implies that the penetration is not 
a function of either impact energy, ~ p b h, or impact mo- 
mentum, ~ ph/i 1 / 2 . 

The penetration of projectiles into granular media has 
also been measured recently by other groups 0, m> ■ 
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Ciamarra et al. |l7j performed quasi-two dimensional 
experiments in which a steel cylinder was dropped side- 
ways into a packing of rods. The impact speeds varied 
by about a factor of five, and the penetration depths 
varied from about 1.5 to 7 times the projectile diame- 
ter. They report that the projectile deceleration is time- 
independent and proportional to the impact speed. This 
implies that the stopping time is constant and that the 
penetration depth isproportional to the impact speed, 
de Bruyn and Walsh |18j performed experiments in which 
two different diameter steel spheres were dropped into 
glass spheres of five different bead sizes. The impact 
speeds varied by about a factor of four, and the penetra- 
tion depths varied from about 1.2 to 5 times the projectile 
diameter. They report that the penetration depth is lin- 
ear in impact speed, but with an intercept d a that can be 
positive or negative. This is modelled in terms of a Bing- 
ham fluid, where the granular medium exerts a force on 
the projectile according to a yield stress and an effective 
viscosity, F(v) = —F — bv. Negative intercepts for depth 
vs speed arc predicted by this model. Even more recently 
Lohse et al. |l9j | performed experiments in which spheres 
are dropped at zero free-fall height (h = 0), just barely 
touching the sand. The projectile densities varied widely, 
at fixed diameter, giving penetration depths d a ranging 
from about 1/4 to 6 times the ball diameter. They re- 
port that the minimum penetration is linear in projectile 
density. This is modelled in terms of Coulomb friction, 
where the medium exerts a force on the projectile pro- 
portional to its depth, F(z) — —kz. Including gravity, 
this law predicts the penetration depth for nonzero drop 
heights to be d/d = (H/d ) 1/2 . 

To summarize published results for the dependence on 
drop height, d ~ i? 1 / 3 was reported in our first paper 
P|, whereas d ~ v ~ ft, 1 / 2 was found in Ref. E3 an d 
d — d ~ h l l 2 was found in Ref. E3- By comparison 
with the projectile diameter, the penetrations are shal- 
low in Ref. 0j but deep in Refs. E3>E3] Thus there may 
be no conflict; the experiments could simply fall into dif- 
ferent scaling regimes. However, evidence of d ~ H 1 ^ 3 for 
deeper penetrations of cylinders is reported in 5]. Fur- 
thermore, Ref. 18] states that our shallow penetration 
depth data of are well-described by their model. This 
raises the possibility that our respective data sets are ac- 
tually consistent, and that one of us is mistaken as to the 
specific power-law behavior. 

The published results for the dependence on projectile 
density are also not in agreement. Our scaling law im- 
plies d ~ p& 3 / 4 for h = 0, Eq. ©, whereas d a ~ pb is 
reported in Ref. |l9j|. Again the penetrations are more 
shallow in our work, so the respective experiments may 
simply be in different scaling regimes. Furthermore, our 
beads are large enough that grain-grain cohesion is neg- 
ligible, whereas the grain size and packing fraction are 
both considerably smaller in Ref. [ljj. To date, we are 
the only group to report density scaling for h > 0. 

Altogether, the results of Refs. 0, El El El suggest 
that there may be three distinct sets of impact behavior. 



(1) Shallow penetration into non-cohesive media, where 
Eq. J2J holds ; (2) Deep penetration into non-cohesive 
media, where d — d a ~ v a and F(v) = —F — bv hold 
(l7L Hq ; and (3) penetration into small tenuously packed 
grains, where F(z) — ~kz holds [19j. 

In this paper we (1) provide more details for our orig- 
inal Letter |4|, and we (2) report on new experiments 
designed to clarify the experimental situation for shal- 
low impacts. Our approach is both to improve the 
reproducibility and accuracy of the measuring appara- 
tus, and to systematically and widely vary the drop dis- 
tance, the ball diameter, and the ball density. We also 
adopt the gas-fluidization preparation method, to see if 
it changes the scaling. We shall demonstrate that the 
deviations of data from Eq. are mainly statistical. 
We shall also demonstrate that better collapse can be 
achieved by Eq. than by the impact speed scaling 
of Refs. [l7l . Hq |. This reaffirms the correctness of our 
original work, Refs. 0, 0| , and negates the statement in 
Ref. E3 that our results can be described by their model. 



I. METHODS 

Our granular medium is P-0140 A-Series technical 
quality solid glass spheres from Potters Industries Inc. 
(PA). The beads are slightly polydisperse, with a diam- 
eter range of 0.25-0.35 mm as set by US sieve sizes 45- 
60. The quoted density of the glass material is 2.5 g/cc. 
In our previous work 0, , we poured the beads into a 
beaker and then gently swirled and tapped it to achieve a 
horizontal surface and a random close packing fraction of 
about 0.64. In case this led to irreproducibility of pack- 
ing or surface angle, we now prepare the system by air- 
fluidization. The sample container consists of a plexiglass 
tube with 8- inch outer diameter, 1/4- inch wall thickness, 
and 5-ft height. The top is open to air, while the bottom 
consists of a Brass sieve with 90 /im mesh opening. Un- 
der the sieve is a windbox consisting of a plexiglass tube 
of same diameter but 12-inch height. The glass beads 
are poured onto the sieve to a depth of approximately 8 
inches. Dry air is then blown at high rate into the bottom 
of the windbox, and up through the glass beads, until all 
fines and humidity are removed. Prior to each drop, the 
beads are more gently fluidized, and the airflow is grad- 
ually reduced, so that a flat level surface remains. It is 
crucial to turn down the airflow very slowly, in order to 
avoid large gas bubbles that leave behind surface irregu- 
larities. In earlier work |20| we found that this procedure 
gives a packing fraction of 0.590 ± 0.004, as expected for 
hard non-cohesive spheres. 

We employ two series of spherical projectiles. The first 
is wooden spheres of density p\, — 0.7 g/cc and varying di- 
ameter: D b = {l/4, 1/2, 5/8, 7/8, 1, 3/2, 2, 3} inches. 
The second is 1-inch diameter spheres of varying den- 
sity: hollow polypropylene, 0.28 g/cc; wood, 0.7 g/cc; 
nylon, 1.2 g/cc; teflon, 2.2 g/cc; ceramic, 3.8 g/cc; steel, 
7.9 g/cc; tungsten carbide (WC), 15 g/cc. These are held 
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and dropped with zero speed from the center of the sam- 
ple tube using a suction mechanism. In comparison with 
our previous work the new sample dimensions and 
maximum ball diameter are all about twice as great, but 
the maximum drop heights are comparable. Since the 
penetration depth grows less than linearly with ball di- 
ameter, we judge that sample-size effects are negligible. 
Furthermore, we never observe any grain movement at 
the edge of the sample as a result of impact. 

The height of the sand surface, the height of bottom of 
the ball prior to drop, and the height of the top of the ball 
after the drop, are all measured using a micro-telescope 
(Titan Tool Supply, Cathetometer TC-11) mounted to a 
height gauge. The sample container and height gauge are 
approximately 1-foot apart, both resting on an optical 
bench. From the height gauge readings, we deduce the 
free-fall height h, the penetration depth d, and the total 
drop distance H = h + d. This method permits study of 
penetrations no deeper than the ball diameter, since the 
top of the ball must be visible from the side. For slightly 
deeper penetrations, until the ball becomes fully buried, 
we estimate the depth from the height of the suction 
mechanism when it is brought into contact with the top 
of ball. 



II. RESULTS 

Raw data for penetration depth vs free-fall height are 
displayed on a log-log plot in Fig. 2] for three example 
wooden spheres and for all 1-inch spheres. The mini- 
mum penetration depth, d Q for h = 0, where the ball 
bottom was initially just in contact with the sand sur- 
face, is displayed along the left axis. For decreasing h, 
the penetration depths extrapolate smoothly to the h = 
limit. For increasing h, the data appear to approach a 
1/3 power-law, shown in Fig. ^ by solid lines. There is 
no evidence of h 1 / 2 behavior, given by the dashed lines, 
which would correspond to the d ~ v scaling of Ref. • 
Fig. ^demonstrates clearly that d cannot be a power-law 
over the full range of h due to the nonzero intercept, 
do > 0. One possibility for simple scaling is that d is a 
power of total drop distance, H = h + d, as advocated in 
Ref. . Another is that d — d is a power of the free-fall 
distance h, as advocated in Ref. [18| . In the next sections 
we investigate both these possibilities. 



III. TOTAL DROP DISTANCE SCALING 

The total drop distance, H = h + d, is a relevant pa- 
rameter because it relates to the average stopping force 
via Eq. QJ. Thus, in Fig. El we replot all the penetration 
data of Fig.^vs H. Now the minimum penetration data 
points, for free-fall height h — 0, lie along the line d = H; 
no data may lie in the shaded region d > H. Whereas in 
Fig. n the data trended toward d ~ h 1 / 3 for large drop 
heights, now in Fig. |3 all the data lie along d ~ H 1 / 3 




h (cm) 



FIG. 1: Penetration depth d vs free-fall height h for (a) 
wooden spheres of different diameter, and (b) one-inch spheres 
of different density. The granular medium is post-fluidized 
glass beads of diameter range 0.25 — 0.35 mm. The projectile 
materials and densities are as follows: hollow polypropylene, 
0.28 g/cc; wood, 0.7 g/cc; nylon, 1.2 g/cc; teflon, 2.2 g/cc; 
ceramic, 3.8 g/cc; steel, 7.9 g/cc; tungsten carbide (WC), 
15 g/cc. 

power-laws. For the lighter spheres, which never sub- 
merge, the deviation from power-law behavior is purely 
statistical. For these data sets, the dynamic range in H 
is two to three decades, enough to give confidence in form 
and a few percent uncertainty in exponent. For denser 
spheres, the H 1 ^ 3 power-law fits gives an acceptable de- 
scription, but may deviate for penetrations deeper than 
about a ball diameter. 

The power-law behavior of Fig. [3 is further analyzed 
in Fig. |3 where we plot the proportionality constant of 
the power-law fits, d/H 1 ^ 3 , as a function of projectile 
properties. In the top plot, Fig. we display results for 
the wooden spheres vs their diameter. Note that three of 
the points correspond to the three example data sets and 
fits shown in Fig. |2Ji. Evidently, to within statistical un- 
certainty, the penetration depth scales as the 2/3 power 
of projectile diameter with a dynamic range of slightly 
over one decade. The penetration depth thus scales as 
d ~ Df, 2 / 3 H 1 / 3 . This expression is dimensionally cor- 
rect, which suggests that we have empirically uncovered 
most of the physics. In other words, the observed D 2 I 3 
scaling lends support to our claim of H 1 / 3 scaling. 

In Fig. |2t> we display the proportionality constant of 
the power law fit, d/H 1 ^ 3 , for all the one-inch spheres 
as a function of their density. Each point corresponds to 
one data set and fit in Fig. Evidently, to within sta- 
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FIG. 2: Penetration depth d vs total drop distance H for 
(a) wooden spheres of different diameter, and (b) one-inch 
spheres of different density. The symbol code is the same as 
in Figure 1. The solid lines are the best fits to d oc H 1 ^ 3 . 
Note that the shaded region, d > H , is forbidden. 



tistical uncertainty, the penetration depth scales as the 
square-root of projectile density with a dynamic range of 
over one and one half decades. We find the same density 
dependence as the jet penetration formula |lfj| . 

Before closing this section, we offer an alternative 
means of analyzing penetration data in terms of total 
drop distance. As noted in the introduction, Eq. J5J) can 
be recast as d/d = (H/do) 1 ^ 3 , where H = h + d, h is 
the free-fall height, and d is the minimum penetration 
depth for h = 0. Thus in Fig.^we check for data collapse 
by plotting d/d vs H/d 0l using measured values of d . 
The scatter of data is not negligible, but the average is 
well-described by the expected 1/3 power-law, shown as 
a solid blue curve. Even tighter collapse onto this curve 
can be achieved if d is treated as an adjustable param- 
eter. For comparison, the d/d = (H/do) 1 ^ 2 power-law 
predicted by the model of Ref. E3 is shown as a dashed 
green curve. For both this model and our observations, 
d is the crucial length scale characteristic of a particu- 
lar system of projectile and granular medium. The value 
of d is proportional to the projectile diameter, Db, and 
a power of the projectile:medium density ratio, Pb/p g - 
The inset of Fig. 0] shows our data for d vs pb, for all 
Df, = 2.54 cm spheres. The results are consistent with 
our expectation, d ~ p^ 4 Eq. ©. For comparison, the 
d ~ Pb observation of Ref. [l!j is shown as a solid green 
curve. We speculate that the small particle size and the 
very tenuous packing in Ref. |19| are responsible for the 
different behavior. 



(a) wood 
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FIG. 3: Scaled penetration depth vs (a) projectile diameter 
and (b) projectile density. Each point corresponds to a fit to 
d oc H 1 / 3 as shown in Fig. |2] The solid lines are the best 
power-law fits, giving exponents of (a) 2/3 for diameter and 
(b) 1/2 for density. The inset of (b) shows the minimum 
penetration depth d a , for h — and Df, = 2.54 cm, as a 
function of ball density along with power-law expectations 
based on Eq. (solid blue line) and on Ref. [l9|| (dashed 
green line). 
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FIG. 4: All the penetration depth data of Fig. Q vs total 
drop distance, both scaled by minimum penetration depth. 
The data all collapse onto d/d a = (H/do) 1 ^ 3 (solid blue 
curve), expected from Eqs. (1-2). The characteristic length 
scale, do, scales with projectile diameter and the 3/4 power 
of projectile density (inset, solid blue curve). By comparison, 
d/do = (H/d ) 1/2 and d a ~ Pb are predicted by the model of 
Ref. EH (dashed green curves). 
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Altogether Figs. 12141 show quite convincingly that the 
penetration depth scales as 



Pb 



(4) 



in accord with Eq. (0). The demonstration here is 
stronger than in our prior work, Ref. Q , because the dy- 
namic ranges are larger and the statistical uncertainties 
are smaller. But more importantly, the demonstration is 
stronger than in Ref. Q because here the projectile diam- 
eters and densities are varied more systematically, with 
one held fixed while the other is changed. Nonetheless, 
Ref. still complements the present work in that it es- 
tablished the dependence of penetration on the properties 
of the granular medium via Eq. © as d~ l/(nPg 1/2 )- 



IV. IMPACT SPEED SCALING 

In Refs. [rH ITsf. penetration depth data are reported 
to scale according to impact speed rather than by Eq. J2J. 
Furthermore, Ref. [18| reports that our earlier data also 
can be scaled by impact speed. Therefore in this section 
we attempt to analyze our new data similarly. We begin 
with Fig. [3 where the penetration depth data of Fig. 2] 
are replotted as a function of impact speed, v — \J2gh. 
Contrary to the suggestions of Refs. our data do 

not lie along the best line fits to d = d ov +v t. For lighter 
spheres the data all curve downwards, while for denser 
spheres the data all curve upwards. Over a subset of 
speeds, e.g. 50 cm/s < v < 400 cm/s as in Ref. 0, the 
fit to a straight line is satisfactory to within experimental 
uncertainty. 

To see if we can make sense of the displayed fits to 
d = d ov + v t 1 we examine one of the fitting parameters 
in Fig. as a function of projectile properties. The top 
plot in shows (d — d ov )/h 1 /' 2 , i.e. the fitting parameter r 
times \/2g, as a function of ball diameter. A reasonable 
power-law fit can be made to D^ 2 ^ , which curiously is 
the same diameter exponent as Eq. (J2J). Though this fit 
gives a fine description, it is not dimensionally simple. 
If (d — d ov ) <x D 2 ^h l l 2 is true, then there must be 
another important length scale in the problem that enters 
the proportionality constant. It would have been simpler 
had we found (d— d ov ) jh 1 ! 2 cx Di, 1 ^ 2 , shown by a dashed 
line. Such behavior clearly differs from the data. 

The bottom plot of Fig. shows the fitting parameter 
t = (d~d ov )/v as a function of ball density. The results 
curve downwards, and cannot be very well described by a 
power-law. Nonetheless, the best power-law fit would be 
to pb 1 / 2 , as shown by the solid line. Curiously, the density 
exponent is the same as in Eq. @. If true, this corre- 
sponds to scaling with the square root of impact energy, 
(d — d ov ) ~ \J pbv 2 . Scaling with impact momentum, 
suggested in the abstract of Ref.0], would correspond 
to (d — d ov ) ~ p b v as shown by the dashed line. Such be- 
havior is vastly different from the data. Contrary to the 
abstract, however, the final scaling advocated in Ref. 0] 
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FIG. 5: Penetration depth vs impact speed for (a) wooden 
spheres of different diameter, and (b) one-inch spheres of dif- 
ferent density. The symbol code is the same as in Figure 1. 
The solid lines are the best fits to d — dov + v t where both 
dov and r are fitting parameters. 
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FIG. 6: Scaled penetration depth vs (a) projectile diameter 
and (b) projectile density. Each point corresponds to a fit to 
d = dov + v T as shown in Fig. |^| The solid lines are the best 
power-law fits, giving exponents of (a) 2/3 for diameter and 
(b) 1 /2 for density. The dashed lines represent dimensionally- 
simpler expectations, but are not consistent with the data. 



6 



is (d— d ov ) ~ Pb~l 2 v . This corresponds to the solid curve 
in Fig. |Ht>, which still is not a satisfactory fit. Even bet- 
ter powerlaw fits can be made in both Figs. Et-b if the 
last point is omitted; however, the resulting exponents do 
not lead to dimensionally simple scaling. While Figs. |SJ- 
El alone do not unequivocally rule out scaling by impact 
speed, the contrast with scaling by total drop distance in 
Figs. 12131 is striking. 

Altogether, the best description of our new data in 
terms of impact speed is 



d-d^ocp^D^h 1 / 2 , 



(5) 



where the intercept, d ov > 0, is a free fitting parame- 
ter as yet unaccounted for. According to the model of 
Ref. [18j , the intercept can be explained by a yield stress 
but only if it is negative, which is not the case for our 
experiments. If Eq. (jSJ is true, the combined density and 
free-fall height dependencies imply that impact energy is 
the crucial parameter, not the impact momentum. For a 
complete understanding, one would still have to account 
for both the free parameter d ov as well as an additional 
length scale in the proportionality constant. 



V. DATA COLLAPSE 

As an alternative means to compare the relative qual- 
ity of candidate scaling descriptions, we now attempt to 
collapse the penetration depth data of Fig. ^ according 
to both the total drop distance H as well as according to 
the impact speed v — \/2gh. For collapse via Eq. 
we plot penetration depths vs ' pgY^D^^H 1 / 3 in 
Fig. Ek- Though the grain density has not been explic- 
itly varied here, we assume the same dependence as ob- 
served in Ref. [^|; this renders the a;-axis dimensionally 
correct. Evidently, in Fig.[7K, the data collapse tightly to 
a straight line with a statistical deviation that is roughly 
a constant percentage. Assuming a proportionality con- 
stant of O.lA/p, as in Eq. J5J), we find that the repose 
angle of the post-fluidized glass beads is 9 r — 21°. This 
is slightly smaller than the value measured in Ref. |4( 
for glass beads at random close packing, 6 r = 24°, as 
expected. 

For collapse via Eq. J3J), we subtract the fitted 
intercept from the penetration depths and plot vs 
(p b /p g ) 1/2 D b 2/3 h^ 2 in Fig. 03. The degree of collapse 
is noticeably not as tight as in Fig. [7^, with a percent- 
age deviation that blooms for smaller penetrations. Re- 
latedly, the systematic curvature of the data away from 
(d — d ov ) cc ft, 1 / 2 , seen in FigOJ), is reflected here by 
a deviation of the actual penetration depth for h = 
from the value of the fitting parameter d ov . If instead 
we plot d — d , where d is the actual observed pene- 
tration for h = 0, then the degree of collapse in Fig. [7|3 
is notably worse. Also, consistent with the diameter and 
density dependence shown in Fig. El the degree of collapse 
is worse in Fig. when the penetration depth is plot- 
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FIG. 7: Penetration depth vs (a) scaled total drop distance 
and (b-d) scaled free-fall height. In (a), the exponents for 
sphere density, sphere diameter, and total drop distance are 
taken from the fits of Figs. 121-11 In (b), the exponents are 
taken from the fits of Figs. 15161 In (c), the exponents are 
taken according to the final prescription of Ref. [lS( |. In (d), 
the exponents are taken such that the x-axis is a length scale 
proportional to impact momentum. In (b-d), d ov is a fitting 
parameter from Fig. that is systematically different from 
the penetration depth at zero impact speed. 



ted vs (pbl Pg) 1 ^ 2 D^^h 1 / 2 , the final scaling advocated 
in Ref. Eg. 

Finally, in Fig. |7|i, we make one last attempt at col- 
lapsing our data. The abstract of Ref. 18] states that 
"...the penetration depth of the spheres increases linearly 
with the incident momentum of the projectile, but with 
a zero-momentum intercept that can be positive or neg- 
ative." According to this prescription, we subtract the 
fitted intercept and plot the otherwise-raw penetration 
depth data of Fig.^vs the dimensionally-simplest quan- 
tity proportional to momentum: (pi,/ p g )Dt, 1 / 2 h 1 / 2 . This 
gives a nearly-random scattering of data points without 
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the least hint collapse. Thus, impact momentum does 
not determine the penetration depth for our data. 

VI. CONCLUSION 

Our new data for the shallow penetration of spheres 
into a loose granular medium strongly support our previ- 
ous conclusions, Eq. for the scaling of penetration 
depth. By improving preparation reproducibility and 
measurement accuracy, we demonstrate that the depth 
scales as the 1/3 power of the total drop distance H. 
In particular there is a nonzero penetration even for zero 
drop height h, where the penetration depth equals the to- 
tal drop distance. By systematically varying the projec- 
tile diameter at fixed density, and by systematically vary- 
ing the projectile density at fixed diameter, we demon- 
strate that the depth scales as the 1/2 power of pro- 
jectile density and the 2/3 power of projectile diameter. 
And by changing the sample preparation from random 
close packing in Ref. Q to a random loose packing here, 



we demonstrate that sample preparation plays no crucial 
role. As long as the medium is loose and noncohesive, 
Eq. J2J) should apply though with a value of n that re- 
flects the packing state. The burning question is now the 
nature of the granular mechanics that gives rise to this 
reaffirmed scaling behavior. The force law cannot be as 
suggested in Ref. [Isj . where the impact momentum of 
the projectile and a yield stress for the granular medium 
are crucial inputs. The positive intercept for penetration 
depth vs drop height, and the scaling with total drop 
distance H rather than with free-fall height h, suggest 
instead that the stopping force may depend on the pro- 
jectile's instantaneous depth as well as its speed. Such a 
scenario would be more in line with recent reports of the 
hydrostatic-like nature of the force on an object moving 
horizontally |2l| or vertically 122^ through a granular 
medium. 
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